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In this paper we consider the 14 moments model of Extended Thermodynamics for 
dense gases and macromolecular fluids. Solutions of the restrictions imposed by the 
entropy principle and that of Galilean relativity for such a model have been until now 
obtained in literature only in an approximate manner up to a certain order with respect 
to thermodynamic equilibrium; for more restrictive models they have been obtained up 
to whatever order, but by using Taylor expansions around equilibrium and without 
proving convergence. Here we have found an exact solution without using expansions. 
The idea has been to write firstly a relativistic model, for which it is easy to impose 
the Einsteinian relativity principle, and then taking its non relativistic limit. 

1 Introduction 

The 14 moments model of Extended Thermodynamics for dense gases and macromolecular fluids 
was firstly studied by Kremer in [I], up to second order with respect to equilibrium. The balance 
equations to describe this model are 



where the independent variables are F, Fi, Fij, Fm, Fan and are symmetric tensors. P<:ij>, Pui, 
Pun are productions and they too are symmetric tensors. The fluxes G^i, G^ij, Gkui, Gkuii are 
constitutive functions and are symmetric over all indexes, except for k. 

In ideal gases, we have also the conditions G^i = F^i, Gm = Fm, Gnu = Pun and, moreover, G^ij 

and Gkiii are symmetric over all couples of indexes; then the present case is less restrictive. 

We want that our system ([T]) is a symmetric hyperbolic one, with all consequent nice mathematical 
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dtP + dkPk = , 
dtPiii + dkGkiii = Pill 



dtPi + dkGki — , dtPij + dkGkij — P<ij> , 
dtFiiii + dkGkiiu = Pun , 



(1) 



properties. To this end, we impose the entropy principle; in the next section it will be proved that 
it is equivalent to the following equations 
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where we have taken into account the definition = F^, from which the following compatibility 
conditions holds 
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Moreover, also in the next section, it will be proved that the Galilean relativity principle is equivalent 
to the following two other conditions 

_ dh' dh' f 9h' dh'\ dh' 

° " 'dX^'^'^^''W,^^'^^\dX;:^'-'^''^'^Wj^^^^^^^ 

° - 'dX^' + ^^''dX] + ^^'^ [dKs ''' '' dX~j)^ ''''dX^i ^ 

In other words, we have to found h' and </)'^ satisfying the equations Q and ([1]); after that, ©i-s 
are useful to obtain A, Aj, Xij, Xm, Xmi as functions of our independent variables F, Fi, Fij, Fm, 
Fun. Lastly, ©7-10 will give the constitutive functions Gki, Gkij, Gkuu Gkuii- 
Now, the restrictions ([3]) and dH) have been imposed by Kremer in [1] only up to secondo order with 
respect to equilibrium. In [2j Carrisi and Pennisi have imposed them up to whatever order (for 
the more restrictive model of ideal gases), but by using Taylor' s expansions around equilibrium, 
without worrying about convergence problems. Here we have found the exact solution without using 
expansions. The idea to obtain this result has been the following 

• Firstly we have assumed a relativistic model for which it is easy to impose the Einsteinian 
relativity principle; the results for this model can be found in section 3; 

• After that, we have shown, in sect. 4, how to take the non relativistic limit of the model in 
sect. 3; to this end we have used a methodology which is easy to find for ideal gases because 
in this case we have suggestions from the kinetic theory of gases; we have adopted this 
methodology also for our more general case. Obvious, the validity of this assumption has to 
be tested at the end by verifying that our results satisfy truly the eqs. ([3]) and (|4]). We will 
be obtain that eqs. ^ are identically satisfied. Instead of this, the condition ([3]) will have 
still to be imposed; it is the only condition for which there is no correspondence between the 
relativistic case and the classical one. 

• Then, in sect. 5 we have taken effectively this limit and found that h' and (p'f^ are determined 
by the following eqs. ([5]) and ([6]), except for 4 scalar arbitrary functions Hq, Hi, H2, H3 
which depend on the scalars 



n 



There remain the further condition the requirement of convexity for the function h' and the 
problem of subsystems. We have exploited them, but don' report here the results, for the sake of 
brevity. We assure only that we have found the exact solution of ^ without using expansions. 
We close this section by reporting the results of sect. 5; in this way it will be not necessary to search 
them throughout the paper and they will be available for the applications. They are 

cp" = hX + HiV^ + H2V^ + H^Vi , (5) 
h' = SHqXi — H1X2 — -H2X3 — , 



with 



1^0' = -2Afc« (6) 
k 4 
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= -2Xl,X,u + -^XuXkaXaU+^XkaXaXpr>ll + 
11 2 

~ l^XiiXkll — XkllXaXall + XkXallXall + 
12 
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— {XaXall)XkbXbll + {trXli,)XkcXcll + -Xll{XaXall)XMl + 
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+ -^XiiX^f^Xhii + -^^XiiXkii — -Xu{trX^^)Xkii + {XabXaiiXbu)Xk + 

2 

— {XabXaXbu)Xkii + -^{trX^f^)Xkii — 2Xi^^Xhii , 

Xi = Xppii , (7) 

X2 = 2XallXall -^XppiiXu , 

= ^XppU {-^X^ - -trA^f, j + 2XabXallXbll — -XuXallXall , 
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X4 = 2Xa^XallXbll - trXabXcllXcll — -XllXabXallXbll + 
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+-^XfiXaiiXaii + SXppu ( --^Xfi + -^XiiitrXli^) - -trXl^, ) , 
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+ -^XiiXallXa — -XiitrX^f, + Y^Xii - 4:XabXaXbll + -tr X^fj , 
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It is interesting that (j)'^ has been determined except for 4 scalar functions Hq, Hi, H2, H^. Instead 
of this, if we have used only the representation theorems without imposing the entropy principle 
and the Galilean invariance, we would have obtained that was depending on 6 arbitrary scalar 
functions, being a linear combination of Xkii, Xkh^hii, ^la-^aii, ^k, ^ka^a, ^la^a- We note that also 
h' is determined in terms of Hq,Hi, H2, H^. These are arbitrary functions of the 8 scalars Xi - X^. 
Instead of this, if we have used only the representation theorems without imposing the entropy 
principle and the Galilean invariance, we would have obtained that all the scalar function are 
arbitrary functions of the following 14 scalars Xu, trXlj^, trXlf^, XaiiXaii, XaiiXa , XaXa, XabXaiiXui, 

Xab^a^blh AafeAaAfc, X'^^XallXbU, X^^XaXbU, A^^AaAfe, XppU, A . 

It is useful now to verify these results: To this end we can substitute eqs. ([5]), ([6]) and ([7|) into ([4]) 
and obtain that they are identically satisfied. The corresponding calculations are long, so it will be 
useful to subdivide them with the following steps. 

• Firstly we can verify that is satisfied with Xi instead of /i', for i = 1, • • ■ ,8; consequently, 
for the theorem on derivation of composite functions, it will be satisfied by whatever function 
of Xi, as h' is. But we have to note that, if we simplify Xg through the Hamilton-Kay ley 
theorem 
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it will become more complicate to verify that is a solution, because it will be necessary also to 
use the identity 

= Sij (^—XahXall^ll + XllXabXallXbll + -XallXalltrXl^ — -XallXallXfj^ + 

+Xij { — XabXallXui + XiiXallXall) — XfjXallXall + XuiXjU -^(^U ~ ^'''^cd) + 

—'^XiiXiniXj)bXbii + XiaXaiiXjbXbii + 2An(iA^-)^AM/ 

which can be easily proved in the reference frame where Xab has the diagonal form. 

• The second step is to verify that ^2 is satisfied in the case Hq = 1, Hi = 0, H2 = 0, = 0. 
Similarly for the case Hq = 0, Hi = 1, H2 = 0, H-^ = 0; then for the case Hq = 0, Hi = 0, 
H2 = I, H3 = and, lastly, for the case Hq = 0, Hi = 0, H2 = 0, H^ = 1. Consequently, it 
will be satisfied for all constant values of Hq, Hi, H2, H3. After that it is satisfied also in the 
general case: In fact, for the property on derivation of a product, the terms where Hq, Hi, 
H2, H3 are not differentiated will simplify, for what above said, and it remains 

]=0 1=1 

for the right hand side of (|3Di written with Xi instead of h'; the result is zero for what already 
verified in the first step. 

2 The principles of entropy and of Galilean relativity 

We want that our system ([T]) is symmetric and hyperbolic, with all the consequent nice mathematical 
properties. To this end we impose that all the solution of eqs. ([1]) satisfy also the entropy inequality 

Ot OXk 

This is equivalent to assume the existence of Lagrange Multipliers A, Aj, Xij, Xm, Xmi such that 

dh = XdF + XidF, + XijdF'^ + XiudF'^^ + XiiudF''^^ , 

d<pk = XdFk + XidGik + XijdGijk + XuidGmk + XuiidGuuk (8) 

besides a residual inequality which we leave for the sake of brevity. The Lagrange Multipliers are 
also called "mean field". Let us now impose the Galilean relativity principle by considering the 
following change of independent variables 

Fiii2...in = ^ ( ^ )"^(ni2---«fe^ifc+i---«n) (9) 
fc=0 ^ ^ 

which can be found in [3] and that, applied to our case, becomes 

F = m, (10) 

Fi = mvi + mi, 

Fij = mviVj + rriij + 2m(jWj) , 

Fill = mm + miiVi + 2miivi + mv^Vi + rmv^ + 2miViVi , 

Fail = rnmi + mv'^ + 4miViv'^ + 2miiv'^ + imiiViVi + 4miiivi. 



Also in [3] we can find how change the constitutive functions Gkii—i„ when the reference frame 
changes, i.e., 



Hkhi^...i„ = ^{A^Kh-ijVi.+^-i^) (11) 
i=o ^-^^ 

where the functions Hki^...i^ are defined by 

Gkh-i„ = VkFh-i„ + Hki^...i„. (12) 

It is interesting that (llip looks like ([9]), except that they don' t act on the index k. In our particular 
case, eqs. (fTT]l become 

Hk = Mk, 

Hki = MkVi + Mki, 

Hkij = MkViVj + 2Mk(iVj) + Mkij , 

Hkiii = MkViV^ + MkiV^ + 2MkiViVi + 2Mkuvi + Mkim + Mkui , 

Hkiiu = MkV^ + AMkiViv"^ + 2Mkiiv'^ + AMkuViVi + AMkmVi + 

+ Mkiiii (13) 

Consequently, the functions Gkii---i„ transform as follows 

Gki = rnviVk + rmvk + MkVi + Mki , (14) 

Gkij = mviVjVk + rriijVk + 2m^iVj)Vk + MkViVj + 2Mk{iVj) + M^j , 

Gkiii = muiVk + miiViVk + 2muViVk + mv^ViVk + nnv^Vk + 

+ 2miViViVk + Mkiv"^ + 2MkiViVi + MkuVi + 2MkiiVi + 

+ MkViV^ + Mkiii, 

Gkiillkk = rriiiiiVk + mv'^Vk + 'irmviV^Vk + 2miiV^Vk + AmuViViVk + 

+ ^mmvivk + MkV^ + iMkiViv"^ + 2Mkiiv'^ + AMkuvm + 



+ AMkiiivi + Mk 



all' 



We note that from (I12p and (llip . for n = 0, and from Gk = Fk it follows Mk = Fk — Fvk- This and 
([I0])i,2 yields Mk = ruk- The new variables m, mj, mjj, mj;/, man and M^, Mjj, M^jj, Mkm, Mknu 
have the same simmetries of Fi^^^^i^ and Gkij^... in- 
Let us now substitute into eqs. ([8]) the expressions which we have above found for the variables and 
the constitutive functions. In this way eqs. dH) become 

dh = X^dm + Af dmi + X^jdruij + Xlndmni + 
+ ^^iiiidmuii + (Af m + 2Afjmj + Xjumadij + 

+ 2Xjiimij + 4Xppggmiii)dvi (15) 

'^('/'fc) = {X^dm + Xldrrii + X^jdmij + Xlndrnm 

+ XiiiidmuiOvk + X^dMk + Xj dMki + X^dMku + 

+ XliidMkiu + XliidMkiiii + (A^m + Afmi + Af^mij 
+ + ^lpqqmiiii)dvk + {X{Mk + 2X-jMkj 

+ XippMkii + 2XippMkii + 4:X^ppggMkiii)dvi, (16) 



with 



A-'^ = A + XiVi + XijViVj + XiUViV^ + XppqqV^ , 

Xf = Xi + 2XijVj + Xiliv"^ + 2XjppViVj + iXppqqViv"^ , 

^ij = ^ij + ^hppVh^ij + 2App(jf j) + AXppqqVjVi + 2Xhhppv'^Sij , 

Xill = Xipp + 4:XhhppVi 

\^ - \ 

^iiii ~ ^ppqq 

The Gahlean relativity principle imposes that h, (pk — hvk, Mf^i, Mf^ij, M^m, Mi don' t depend on 
Vi. For this condition on h and (j)k — hv^ we obtain 



dh 

dvi 



= mAf + 2X(jmj + Xj JmiAj + 2mij) + 
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— = = MkXl + 2MkjXlj+MkuXlpp + 2MHiX\pp + 

+ 4Xl,pggMkiii + h'5ik, (17) 

where h' is defined by 

h' = X^ m + Xlrui + Xlima + Xliirriiii + Xliiimm - h\ 
In this way eqs. (jlSp and (jl6p become respectively 



dh^ = X^dm + Afdmj + Xlidmu + Xlndrnm + Xlmdmiiu (18) 
#i = A-^dMfc + Xl dMki + Af^dMfciz + XiiidMkiii + Xlp^^dUkm 



with (j)^. = — hvk- Let us also define 

= A-^Mfc + XlMki + XliMkii + XiiiMkiii + XlpggMkiiu - 4 ; 

so that eqs. (fTH]) become 

dh' = mdX^ + rtiidXl + mudXli + muidXln + muudXlui 

d(t)'k = MkdX^ + MkidXl + MkiidXli + MkuidXlu + MkmdX^ppqq 



from which, by taking A^, Af , Xjj, X^n, X^-n as independent variables and, by taking the derivatives 
with respect to the various components of the mean field, it follows 

dh' dh' dh' 

dh' dh' 

~ ; 
/ 

(20) 



Mk - ^^^-^^ 

dcp'f^ dcp'j^ 
Mkiii = WTj-j Mkiiii = J ■ 



d>4a dXl, 



These last equations are noting more than ([2]), but in the new reference frame. By substituting 
in (|17p from (j20p we obtain equations whose expression, in the previous reference frame, are 
Consequently, ([2]) and @ are equivalent to the principles of entropy and of Galilean relativity. 
Other interesting aspects can be found in [3], which are adapted for the present case in [5] and 
In order to find the general solution of (j3|), let us write a relativistic counterpart of our equations. 

3 Relativistic extended thermodynamics for dense gases 
and macromolecular fluids 

Let us consider the balance equations 

9^r"^ = , daA'^^^ = I^^ . (21) 

where T"^ isn' t symmetric, while and I^'^ are symmetric only with respect to the indexes 

(3'y. The first of these is the conservation law of momentum-energy, while the trace of the second 
one is the conservation law of mass, so that 

/^^5/37 = 0- (22) 

The entropy principle is expressed in terms of h'^, called (entropy density - entropy flux density) 
tensor, such that 

dah" = a>0 (23) 

for every solution of eqs. (pTI) . 

For Liu Theorem p'J eq. (|23p is equivalent to assuming the existence of the Lagrange multipliers 
A/3, Xjs-y such that 

-a- Xpd^T^P - A/3^(a„^"'3^ - = (24) 
for every value of the independent variables. By differentiating it becomes 

d/i" = A/jdT"^ + Xp^dA""^^ ■ -a + Xp^I^^ = 0. (25) 

If we define h"^ by 

/i° = + XpT'^^ + A/3^A"^^ (26) 

then eq. (f25|l can be rewritten as 

dh"" = T'^^dXp + A'^^^dXp-^. (27) 
This last equation, by taking Xp, Xp^ as independent variables, becomes 

T-' = ^ , ^"'^^ = 1^ (28) 

It follows that from the knowledge of /i'" we obtain T"^ and A^l^"' ; we have only a condition on /i'" : 
it has to satisfy the Einsteinian relativity principle. For well known representation theorems as [8] 
and [9], we have that 

/i'" = /ioA" + /iiA^^A^, + /i2 A "^A^ + /i3A"^A^ (29) 



where the following definitions have been used 



with hi scalar functions depending on 

Q^~^fi ' ^2 =A '^^ ga-i , Qs =A ""'gay , Qa =A °'"' gay , (30) 

Po = A/3A^ , Pi = A^A^A^^ , P2 = A;3A-,A^^ , P3 = A/3A-, A . 

4 The non relativistic limit of the previous model 

In order to take this limit, let us consider a modified procedure of that used in [10], [11] for ideal 
gases. We assume that this procedure holds also for macromolecular gases; this assumption doesn't 
lead to wrong results because we have already verified, at the end of sect. 1, that these results are 
correct. Let us subdivide the procedure in two parts. 

4.1 A first transformation in 3-dimensional form 

Eq. (f2T]l i for (3 = 0,i becomes 

-dtT^° + SfcT'^o = , (31) 
c 

^dtT'^' + dkT''' = . (32) 



Similarly, eq. ([2T]) 9 for (3^ = 00, (3^ = Oi, /?7 = ij becomes 



c c 
Now we adopt the following change of variables 



c 



(33) 



TOO = m^cFa, T'^o = m^G^, T^^ = m^Fl T^' = "^Cf (34) 



^000 ^ mgc^Fg, ^ooi ^ mlcFl A^^i = mlFl^ (35) 

A*^oo ^ ^5^G^^ A^^' = mlGf, A^'^ = ml^Gl'^ (36) 

/°o = cPm^, = g*m^, = ^p'^ml . (37) 

In their terms eqs. (f3T|) . ([32]) and ([33]) become 

9tF2 + 5fcG^ = , 9tF2^ + a^G^* = , OtF^ + dkCl = R (38) 



n 



The last of these can be also subdivided in 
while eq. ([23]) can be subdivided in 

-dth"^ + dkh^ = a . 
c 

With the following changes of names 

/lO = mlh, = a = 

c c 



eq. (poj) becomes 

'm -h Ok 



dth + da'' = a* . 



4.2 A first transformation of the Lagrange Multipliers 

EQ. (p5|) i for a = becomes 

dh^ = Aodr°o + AidT* + Aood^°°° + 2Aoid^°°' + Ayd^o^^' . 

By using eq. ([^Tjl i and P^ - ([55|) . it transforms into 

d(mg/i) = Xod{m^cF2) + Aid(m^F2*) + Aood(m-oC^^3) + 
+ 2Xoid{ml,cF^) + A,jd(mgF3*^') 



from which 



dh = Xomocd{F2) + Ximod{F^) + Xoomlc^diF^) + 
+ 2Xoimlcd{F^) + Xijmld{F'^^) = 
= idF2 + lidFi + r/dFa + /u^dFj + /^ydFg^^' , 



I = XorriQC, £i = Ajmo , rj = Xoouilc^ , m = 2Xoimlc, ixij = Xij-ml . 



with 

Similarly, eq. ([25]) i for a = /c becomes 

dh'' = X^dT''^ + X^dT''' + Aood^''°° + 2Aoid^^°^ + X^jdA'''^ . 
By using eq. ([^T|) and ([M]) . (p6|) . it transforms into 



+ 2Xoid{mlGf) + A,,(i(m5iG^'^) 



which can be rewritten also as 



dcj)'' = Xocmod{G2) + XimQd{Gf) + Aooc^?Tiod(G|) + 
+ 2XQ,cmld{Gf) + Xijmld{Gf^) = 
= IdG'^ + lidGf + r/dGg + ^i.dGf + ^l^jdG'f^ , 



1 n 



where we have used (j46p . The equations (j45p and (j48p represent the entropy principle for the system 
psp : moreover, eq. (j46p gives the transformation of the Lagrange multiphers. 
Let us now find the transformation of h'" ; eq. (|26p for a = is 

/jO = -/I'O + Aor°° + AiT^* + AooA°°° + 2Aio^°*° + A^^ . (49) 

By using eqs. ([MD and dH])-®, eq. ([MD becomes 

mlh = + £mlF2 + ^im|]Fj + r^m^Fg + /iimgF| + /iij-mj^Fg'-^'. (50) 

This can be rewritten as 

h = ^+m+ dFi + vFs + fiiF^ + fiijF;,' (51) 
rriQ 

or 

h = -h' + iF2 + i^F^ + r/Fg + fiiF^ + /x,jF3^^' (52) 

where we have defined 

/I'O = m^h' . (53) 

In this way we have found the counterpart of (I4ip i for /i'". Let us now find the transformation of 
h'" for a = k; from eq. (j26p we find 

/i'^ = -/i'^ + AqT'^o + AiT^^ + Aoo^^°° + 2AioA'''° + A,^^*^*^' (54) 

This, by using again eqs. (f4T]) 9. (l46]) and p4l) - (i36l) . can be written as 

c rriQC mo c 

m^c^ ^ ^ m^c c ^ 

or 

with 

h"^ = Hi^'k^ (56) 
c 

Eqs. (f53l) and (f56l) are the counterparts of (f^TIl i 9 for /i"^ and /i''^. Eqs. (f52]l and ([55]) are the 
counterpart of eq. (j26p for the system (j38p . Let us finish by considering the counterpart of mass 
conservation ([22]), that is + P-'iJjj = or, by use of ([37]), - cRml + ^p'^m^dij = 0. In other 
words, 

4.3 The second transformation in 3-dimensional form 
4.3.1 Suggestions from the kinetic theory for ideal gases 

For ideal gases the variables ^2,^3,^3,^3, have counterparts in statistical mechanics where they 
are defined as moments of the distribution function /, i.e., by means of the following integrals 



F2 = J h^du, Fi = Jh'u'dn, (58) 
F3 = JfYdu, Fi = Jfyu'du, F^^ = J fYu^u^du. 



where 7 is the Lorentz factor (1 — ^)~2. 

If we take the hmits of these expressions for c — > cx) we obtain F2 = F^, F2 = F^ so that they 
will be no more independent variables; to avoid this problem, we take suitable invertible linear 
combinations of the equations before taking the limits. 

In particular, as first equation we take the linear combination of ([Ml) 3 and (j39P i trough the coeffi- 
cients 1 and respectively; so we obtain the following eq. (j69|) i which is the conservation law of 
mass with 

F = F3-^F^ G'' = Gl-\Gf- (59) 



moreover, we have taken into account eq. (j57p . 

As second equation we take (ISSP o which can be written as the following eq. (1691) 9 which is the 
conservation law of momentum with 

F' = Fi , G^' = Gf ; (60) 

As third equation we take the linear combination of (f38]) i . ([SHJs, (f39l) i trough the coefficients 2c^, 
— 2c^, 2 respectively; so it becomes 

+ dkG^^^ = (61) 

which is the conservation law of energy, with 

= 2c2(F2 -F3 + ^Fi'), G^^" = 2c2(G^ - + ^Gf ) . (62) 



As fourth equation we take (I39p 9 which we write as 

dtF<''^> + dkG^<'^> = p<'^> (63) 

where 

F<'^> = F<'^> , G^<'^> = gI<'^> . (64) 

We transform furtherly eq. (|63p adding to it eq. (|6ip multiplied by ^ and obtaining the following 
eq. (I69])3 with 

F'' = Fl' + ^(2i^2 - 2F3 + ^) (65) 
G^'^ =gI'' + ^{2G^2-2G'l + ^) (66) 



Equation (j69p q encloses both (j6ip and (j63p . As other equation we take the linear combination of 
psp ^ and ()38p 9 trough the coefficients 2c^ and — 2c'^, respectively; so we obtain the following eq. 
([69]) /! with 

pill ^ 2c^{fI - Fi) , G^*" = 2c2(G^' - Gf) , /' = 2c^Q^ . (67) 

Finally, as last equation we take the linear combination of (j38P i . (j38p -^. (j39P i trough the coefficients 
— Sc"^, +8c^, — 4c^ respectively; so we obtain the following eq. (169]) ^ with 



pan ^ _8c4i?2 + Sc'^Fs - ic^F^^ 
Qkiiii ^ _Sc'^Gfl + 8c^G§ - Ac'^Gf (68) 

piiii ^ 8^4^ _ 4^2^// ^ 8g2p«« _ 4^2^// ^ 4^2^^ ^ 



1 o 



where eq. (jSTj) has been used. The complete system is 

dtF + dkG^ = , dtF' + dkG^' = , dtF'^ + dkC'"'^ = p<'^> , (69) 

Note that the only difference between this system and ([T]) is that in (|69|) i intervenes G'^, while in 
([T])i there is F^'; for this reason we will have to impose the further condition = F'^ . 

4.3.2 Reasons for the above choice of coefficients. 

The coefficients in the above linear combinations have been chosen for the following reasons: from 

dSH]), (I59l)i, iQDi, (iMDi, dell, dSIDi, (EHDi it follows 



F = j h\l -^)d , r = l h'u'duu = I h'du , = 2c' I h'{^ - 1 + ^)du 
F<'J> = J h\u'u^ - ^u'5'^)dvi , F*" = J h'^u'il - ^)du , 

J 7 

which have limits J fdw, J fu^du, J fu'du, J fu'^^u^^du, J fu^u'du, J fu^dn, where we have taken 
into account also that 

Prom the previous expressions it follows that F^^ has limit J fu^u^du. Obviously, the previous 
properties hold in the case of ideal gases; we assume that the corresponding change of equations is 
appropriate also for dense gases and for macromolecular fluids. In this way generality is not lost 
because the Galilean relativity principle has the same form for both cases. 

4.4 Second Transformation of the Lagrange multipliers 

The change of equations in the previous section induces another one on the Lagrange multipliers 
and we want now to determine it. 

We observe that the equations (lMl)i,(l6Q])i, ([65]), (l671)i, ([681)1 give F, F\ F'^ , F*", F"" in terms of 
F2, F|, F3, F3, Fg-'. Let us now take the inverse of these relations, which are 

P2 = ^ + F , F^ = F' (71) 

pppll pll pill . . sjij pppll 



and these we now substitute in eq. (jl5|) : so it becomes 



pll pppll pll 

dh = ld{^+F)+kd{F')+r,d{F + ^ + -^) + 

pill 1 pppll 

+ f,,d{— + r) + fi,,d{F^^ + -5'^ 



XdF + XidF' + X.jdF'^ + Xiudr''XppudFPP'' (72) 



with 
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X = l + r] , Xi = li + fii , \j = IJ-ij + {-^ + , (73) 

\ii - TTo > ^ppu -t:i + 



2c2 ' 4c4 ' 12c2 ■ 

In this way we have found the first part of the entropy principle for the new system. 
For the sequel it wih be useful to take the inverse of eqs. (j73|) . They are 

■q = Sc^Xppii - ^c^Aiz + A , i = -Sc^Xppii + ^c^^ii , /"j = 2c^Aiii , (74) 

h = K — '^(?\u , = Kj — (4c^App;i ^ ~'~ c^'^'^*-' ' 

Similarly, eqs. §9^2, &2, §8^2 gi^e G'', Qkij^QkUi^ Qkuii ^ -^^ ^^^.^g Qk^ Qki^ 

G|, G3* G^-' . The inverses of these relations are 

r~<kU 

G^ = — + G^ , G^' = G^' (75) 

f-ikppll r^kll /~ikill -i /^kppll 

Gl = G^ + '^ + ^ , Gf='^ + G^^ , Gf = G^^^ + -5^=^ , 
which now we substitute into eq. (j48p which now becomes 

/-~ffc// rtkppll rtkll 

dcj)^ = idi^ + G'')+£^dG''' + 7]d{G'' + ^^^ + ^)+ (76) 

4c^ 



f-ikill 1 r~<kppU 



= XdG^ + XidG''' + XijdG'''^ + A.ijdG'^*" + Appi^dG'^PP" , 

thanks to eq. ([73|) . 

In this way we have obtained the second part of the entropy principle. We deduce now the trans- 
formation of h' and 0"^; to this end, let us take eq. ([52]) and substitute eqs. ([71]) : so we obtain 

h = -h' + li^ + F)+kF^ + r^{F + ^ + ^)+ (77) 

pill xij pppll 

which, for eq. (|73p . can be written in the following way 

h = -h' + XF + XiF' + XijF'^ + Xiiir'' + XppuFPP'' (78) 
From eqs. ([55]) . by substituting from eqs. ([75]) we find 



/~<kll 

b"' + + G^) + iiG''' + (79) 

(^kppll (J^^^ Qkill ^ij Qkppll 



4c4 ' c2 ^ ' '^'^ 2c2 ' ^ ' ' ^'^'^ ' 12c2 ^ 
= + AG'^ + AiG^^ + XijG^'^ + XiiiG^^^^ + XppuG^PP^^ , 

for eqs. ([73]) . 

The equations ([72]) and ([75]) gives the entropy principle for the balance equations ([M]) . The equations 
([78]) and ([79]) give the counterparts of ([52]) and ([55]) for the new system ([69]) . We have now to obtain 
h' and (j)'^ from ([29p expressing them in terms of the new Lagrange multipliers and then taking the 
limits for c — > 00. The functions h' and (j)'^ are called "potentials". 
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5 Determination of the potentials 



With eqs. (j78p and (j79p we have found that the functions h' and (j)'^ for the new system (j69p are 
exactly the same obtained, through the above mentioned passages, from the 4-potential h"^ for the 
system (j2ip . Consequently, we may obtain them from the result (j29p of the relativistic system. To 
this end it is necessary to write the relativistic Lagrange multipliers in terms of those for the system 
([69]) . By deducing Aq, Aj, Aqo, Aqi and Ajj from (06]), and by substituting in their expressions eqs. 
([Tip we obtain 



- — 

7 2 



A'' 



+ 



c 

mo 




+ - 

-4:\ppllSij/ \Kll Oi 



1 il^ii 



-A 0, 



Oj Ajj + \\u5ij 



+ 




+ 



1 



-2A,7, 



Oi 



+ 



(80) 



We can now begin to evaluate the 4- vectors intervening in (j29p and the scalars (j30p . 



5.1 The scalars Qi - Q4 . 



Let us begin with (|30P i . From eq. (jSUP i we find 



Ql = — ( -20App;; + -^-^kl 



(81) 



Consequently, if we assume that the scalar functions depend on Qi as composite functions through 

2 

Y^Qi, then their limits will be functions of Xppu = Xi. 



-20c 



Let us consider now eq. (j30|) 9. From eq. (jSUP i we find 



Q 



4r 



112 a: 



+ 0(1; 

c 



where O(-) has limit zero when c goes to infinity. Therefore, if we assume that the scalar functions 

4 

depend on Q2 as composite functions through ^Q2, then their limits will be functions of 112A: 



ppii ■ 



But this is not independent from Xi, so that this result is too much restrictive. The idea is now to 

4 

find a number k such that ^{Q2 + kQl) has zero limit for c going to infinity. We find k = 
After that, we see that 



16 1 

-'^Ka^aii + —Xppiiki + O(-) 
5 c 



(82) 



Therefore, if we assume that the scalar functions depend on Qi and Q2 as composite functions 

2 4 

through Ql and —t{Q2 — mQi)^ then their limits will be functions of Xi ed X2. 



1 c 



Going on in a similar way, we look for the numbers ki and k2 such that ^[Qs + kiQi{Q2 — ^Qi) + 



k2Qi] has zero limit for c going to infinity. We find ki 



Q3 - rrTQi{Q2 - l^Qi) - -r^Ql = — 
1(J 25 125 ttIq 



— jq and k2 
2 / 3 



25 

j^. After that, we obtain 



5^3 + 0(- 



(83) 



Then the limit of a scalar function depends on Xi, X2 and X3. 

Similarly, we search the numbers k^, k4, k^, k^ such that ^{Qa + ^sQiiQa — jqQi{Q2 — ^Qi) ~ 
mQi\ + ^4(Q2 - + hQl{Q2 - ^Qi) + hQf] has zero limit for c going to infinity. We find 



16 
15' 



2' 



Q4 



16 

15 



Qi 



14 
'25' 



Q-s 



19 
625- 



After that, we obtain 



14 

25 



25 



10 



Ql 



Qi 



Q2 

19 
625 



25^^ 



11 
125 



Q2 



-Ql 

25^^ 



(84) 



-2X4 + 0(-) 
c 



Therefore, the limit of a scalar function depends on Xi, X2, X3 and X4. 

Before proceeding with the other scalars ([SDJs-g, let us evaluate the 4- vectors of which 

linear combination. 



IS a 



5.2 Limits of the 4- vectors in the expression of the relativistic 
potential 



From the eqs. (|53|) . (j56|) and (j80j) 9 it follows that the term /iqA" contributes to (p the term —2HQXkii 

3 

and to /i' the term SHoXppu, where Hq is the limit of -^/io for c going to infinity. 



But, from eqs. (1561) and (I80p it follows that the term /iiA^'^A^ contributes to the term -^/ii 16A„p;;Afci; 
which is parallel to the previous one. In order not to lose generality, it is better to look for a number 
a such that (A"'''A^ + aQi A°) ^ has zero limit. We find a = — |. After that, we obtain 



QiA^ 



m: 








(85) 



2 3 

Now a linear combination, with arbitrary coefficients, of A", A^'^A-y, A "'^A-y, A^'^A-y is also a linear 

2 3 

combination, with arbitrary coefficients, of A", A"'^A^ — 5 Qi A", A "'*^A^, A"'^A^ so we can suppose 
that in ([29]) there is A^'^A ^ 
that the term hi{X"'^X^ — -g 



7 5 Ql A° instead of A^'^A^,; After that, from eqs. (j53l) and (156]) it follows 
^ Ql A") contributes to cp"' the term Hi{-2XihXhii + ^XppuXi + ^XuXm) 

3 

and to h' the terms H1X2, where Hi is the limit of -^hi for c going to infinity. 

rrtQ 

Proceeding furtherly in this way, we search the numbers ai, 02, 03, such that 

^ |a "^A^, + aiQiA^'^A^ + [aaQf + a3{Q2 - ^Qj)] A"), has zero limit for c going to infinity. We 
find oi = 



^' "3 



After that, we obtain 




7 

25 



A° 



(86) 



This 4- vector can replace A "'*^A^ in (j29p : so we find that it, together with the factor /12, contributes 

^kll^a^all + ^kKllKll + 

■/i2 for c going to infinity. 



to (j)' the term H2{—2\^,^\hii + -^XuXka^aii + ^XkaXa\pn — 2^\i^ku 



+(trXlf,)Xkii - -^XppiiXiiXk) and to h' the term H2X3, with H2 hmit of 

u 

Finally, we search the numbers 04, 05, ag, 07, ag, ag, such that 
^ { A "^A^ + a^Qi X "^A^ + [a^Ql + 05(^2 - ^Ql)] A"^A^ 

+ [q^tQs + flsQi + fi9Qi(Q2 — ^Qi)] }i hs-s zero limit for c going to infinity. We find 04 



as = -I, ae = — ^, 07 = — -I, as = og = |. After that we have the result 



4 

'5' 



3 42 

A "^A-, - -Qi X "^X^ + 



-lQ3 + ^Ql + lQiiQ2-^Ql) 




A" 



A"TA^ + 



(87) 



It may replace A ""^A^ in (I29p . so that it, together with the factor h^, contributes to cf)'^ the term 



H3 [2A, 



ppii 



-2Xl^Xh-trXl,Xk 

17 \2 



^XiiXkaXa + ^A^;Afc) + {XkhXh){XauXall) + 



jXii{XallXall)Xk — ^XfiXkaXall — {XallXall) XkbXui + {tr X^i,)XkcXcll + |Az/(AaAai/)Afc// + | A// A|^A/J// + 
| A;; (tr A^ J Afc// + {XabXauXMi)Xk - {XabXaXbu)Xku + |(trA^JAMi - 2Xlf^Xhu] and to h' 



the term H^X4, with limit of -^hs for c going to infinity. 



5.3 The scalars P 

From eq. (f80l) 9 we find Pq 







c6 



^^^ppii + ^(c)j ^^^^ ^'^^ limit of ^-Pq is a function of Xi. In 
order to obtain a less restrictive result, we may substitute Pq with Pq + ^Qi"T'o'^^' because this will 

ppll 



eliminate the term -^64A^„„ di Pq: in fact we obtain 



-2X0 



0(1) 



In this way we have a better result, even if it is still not enough; to this end let us substitute 
Po + ^Qi'itt-qc'^ with Pq + ^QlrriQC^ + 2{Q2 — ^Q\)'ni'^c^ and the result is satisfactory because is 
equal to 



/ 2 

— o I — -A;; + 16App/;A 



m, 







4AaA,z; + 2trA^fe + 0(-) 
c 



In this way we have found the new scalar X5. 
Let us now consider the scalar Pi in ([30]) fi. It is obvious, for eq. 



85|) that it is better to replace it 
with Pi — |<5iPoj but this has limit a scalar which is a function of X\, X2, X^. Briefly, let us look 
for the numbers 61 and 62 such that 



7 



Pi - 7Q1P0 + hiQi ( Q2 - —Qi ) m^c' + 62 



Q3--rQi(g2 



10 



25 



11 
125 



2 2 



1 ^7 



has zero limit. 

We find = 1, ^2 = I- After that we have 



Pi 



IQiPo + ^Qi (q2 - ^Ql) mlc" + ^ 



Q3-^Qi{Q2-^Ql)--^Ql 



7 

25 



11 

125 



2 2 
ITIqC 



= ^(x6 + 0(i: 
mg \ c 

from which the new scalar Xq. 

Let s consider now the scalar P2 in (jSOP y. It is obvious, for eq. ()86p that it is better replace it with 
P2 — \Q\P\ + [^Qi ~ \ {Q2 — ^Qi)] Po'i but also this is not enough; then we search the numbers 
63 and 64 such that 



m, 



+&3Q1 
Q4- 



^ P2 - tQiPi + 



1 



7 



25^1-2 1^^-25 ^1 



Po+ 



16 
15 



Qi 



11 

125 



Ql 



+ 



--^{Q2- l^Ql) - '4qI (q2 -IqD- \ 



25 



19 

625 



has zero limit. 
We find 63 = -4, 64 



1. After that it follows 



P2 - tQiPi + 
5 



25 



4 

+ 15^^ 



Q3 



— Ql {Q2- — 

10^ V 25^^ 



Q2 - T^Qf 

11 



125 



Po + 



2 2 I 



+ Q 



16 
15 



Ql 



14 



10 



-Ql 

25^^ 



11 

125 



Q? 



+ 



from which the new scalar Xj. 

Finally, let us consider the scalar P3 in (|30]l «. It is obvious, for eq. (f87|l that it is better replace 

it with P3 - IQ1P2 + i^Ql - i(Q2 - ^Q?)] Pi + [-IQ3 + ^Qf + |Qi(Q2 - ^Q?)] Po; but this 

also is not enough; then we look for a number 65 such that 



rQiP2 + 



1 



Qi-i;iQ2 



7 
25 



Ql) 



Pi + 



16 



-^Q3 + ^Q? + ^Qi(Q2-^Q?) 



Po\ + 



+65Qi^|Q4--Q 



Qs - ^Qi ( Q2 - ^Ql ) - i^Ql 



7 

25 



11 

125 



+ 



has zero limit. 

We find 65 = i Then we can evaluate 



P3 - TQ1P2 + 
5 



1 



^^1 + 



\qz + ^Q\ + \qi{Q2-^QI) 



P0 + 



16 
15 



Qi 



Q3 



14 



—Qi (q2- — <3? 

10^ V 25^^ 



11 
125 



Ql 



+ 



7 



19 



2„2 



from which the last new scalar Xg. 



NOTE: We have obtained the decomposition (jSOp by considering a modified procedure of that 
introduced in [10], [11] for ideal gases; instead of this, if we use exactly the procedure of [10], [TT] , 
we have to substitute (IHOl) with 



= — 



+ 



A 0,- 




+ 



Another procedure is present in literature ( [12j . |13j ) also for ideal gases. If we want to follow it, 
then we have to substitute (l80l) with 



A^ - — 



A^ 






0, ^ 


1 




^<ij>J 





(89) 



We have performed calculations also with (j88p and with (j89p instead of (j80p . but we don' t report 
them for the sake of brevity. The interesting result is that the pertinent polynomials in 1/c are 
different between them and from those here obtained, but the limits for c going to infinity are the 
same with all 3 approaches! 
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